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ABSTRACT 

The spectrum of Seyfert 1 Galaxies is very similar to that of several Galactic Black Hole Candidates (GBHCs) 
in their hard state, suggestive that both classes of objects have similar physical processes. While it appears that an 
accretion disk corona (ADC) model, where the corona sandwiches the cold accretion disk, is capable of explaining 
the observations of Seyfert galaxies, recent work has shown that this model is problematic for GBHCs. To address 
the differences in the spectra of Seyferts and GBHCs, we consider the structure of the atmosphere of the accretion 
disk in a region near an active magnetic flare (we refer to this region as transition layer). We show that the 
transition layer is subject to a thermal instability, similar in nature to the ionization instability of quasar emission 
line regions. We hnd that due to the much higher ionizing X-ray flux in GBHCs, the only stable solution for the 
upper layer of the accretion disk is that in which it is highly ionized and is at the Compton temperature (kT ^ a 
few keV). Using numerical simulations for a slab geometry ADC, we show that the presence of a transition layer, 
here modeled as being completely ionized, can signihcantly alter the spectrum of escaping radiation for modest 
optical depths. Due to the higher albedo of the disk, the transition layer leads to a reduction in reprocessing 
features, i.e., the iron line and the X-ray reflection hump. In addition, if most of the accretion energy is dissipated 
in the corona, the thermal blackbody component is reduced, giving rise to a lower Compton cooling rate within the 
corona. Although a transition layer does allow for higher coronal temperatures, global two-phase, slab-geometry 
ADC models still cannot have coronal temperatures high enough to explain the data unless the optical depth of 
the transition layer is unphysically large (^ 10). However, a model having a patchy corona rather than a global 
corona appears promising. Thus, it is possible that due to the thermal instability of the surface of the accretion 
disk, which leads to different endpoints for GBHCs and Seyfert galaxies, the X-ray spectra from these two types 
of objects can be explained by a single unifying ADC model. 

Subject headings: accretion disks — black hole physics — Cygnus X-1 — galaxies: Seyfert — magnetic helds 
— radiative transfer 


1. INTRODUCTION 

The X-ray spectra of Seyfert Galaxies and Galactic Black 
Hole Candidates (GBHCs) indicate that the reflection and re¬ 
processing of incident X-rays into lower frequency radiation is 
an ubiquitous and important process. For Seyfert Galaxies, the 
X-ray spectral index hovers near a “canonical value” (~ 0.95; 
Pounds et al. 1990, Nandra & Pounds 1994; Zdziarski et al. 
1996, but also see Brandt & Boiler 1998), after the reflection 
component has been subtracted out of the observed spectrum. It 
is generally believed that the universality of this X-ray spectral 
index may be attributed to the fact that the reprocessing of X- 
rays within the cold accretion disk leads to an electron cooling 
rate that is roughly proportional to the heating rate inside the 
active regions (AR) (Haardt & Maraschi 1991, 1993; Haardt, 
Maraschi & Ghisellini 1994; Svensson 1996). 

Although the X-ray spectra of GBHCs are similar to that 
of Seyfert galaxies, they are considerably harder (most have a 
power-law index of E ~ 0.7), and the reprocessing features are 
less prominent (Zdziarski et al. 1996, Dove et al. 1997). It is 
the relatively hard power law (and therefore the required large 
coronal temperature) and the weak reprocessing/reflection fea¬ 
tures that led Dove et al. (1997, 1998), Gierlinski et al. (1997) 
and Poutanen, Rrolik & Ryde (1997) to conclude that the two- 
phase accretion disk corona (ADC) model, in both patchy and 
slab corona geometry cases, does not apply to Cygnus X-1. 

One of the main problems with the global slab-geometry 
ADC model is that, for a given coronal optical depth, no self- 


consistent coronal temperature is high enough to produce a 
spectrum both as hard and with an exponential cutoff at an en¬ 
ergy as high as that of Cyg X-1 (Dove, Wilms, & Begelman 
1997). However, this result is sensitive to the assumption that 
the accretion disk is relatively cold, such that ~ 90% of the 
reprocessed coronal radiation is re-emitted by the disk as ther¬ 
mal radiation (with a temperature ~ 150 eV). It is this thermal 
radiation that dominates the Compton cooling rate within the 
corona. If the upper layers of the accretion disk were highly 
ionized, creating a “transition layer,” a smaller fraction of the 
incident coronal radiation would be reprocessed into thermal 
radiation (i.e., the albedo of the disk would be increased), and 
therefore the Compton cooling rate in the corona would be re¬ 
duced. Furthermore, as shown by Ross, Fabian & Brandt (1996; 
RFB96 hereafter), the high ionization state of the outer atmo¬ 
sphere of the accretion disk in GBHCs may explain the weak¬ 
ness of observed iron line features in Cyg X-1. 

Recently, Nayakshin & Melia (1997) investigated (via simple 
qualitative considerations) the structure of the X-ray reflecting 
material in AGNs assuming that the ARs are magnetic flares 
above the disk (Haardt et al. 1994; see also Galeev, Rosner & 
Vaiana 1979). They showed that the pressure and energy equi¬ 
librium conditions for the X-ray illuminated upper layer of the 
disk require the gas temperature to be T ~ few x 10® K, leading 
the upper atmosphere of the disk to a low ionization state. 

In this paper, we consider the reprocessing of X-rays from 
ARs for GBHCs. In §2, we show that there should be a thermal 
instability at the surface of the cold disk, causing the temper- 
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ature to climb up to T ^ a few x 10^ K. This temperature is 
roughly the Compton temperature with respect to the coronal 
radiation field. At this temperature the transition layer turns 
out to be almost completely ionized. Recently, Bottcher, Liang, 
and Smith (1998), using an iterative method, a linear Comp- 
tonization algorithm, and the photoionization model XSTAR, 
also found that a highly ionized transition layer should form 
for moderate ionization parameters, i.e., ionization parameters 
thought to be appropriate for GBHCs. In §3, we explore the 
ramifications of this highly ionized transition layer on the ener¬ 
getics of the corona, and investigate how it alters the spectrum 
of escaping radiation. We also discuss whether slab geometry 
ADC models, when transition layers are included, can account 
for the observed spectra of BHCs. In §4, we give our conclu¬ 
sions. 


2. THE FORMATION OF A TRANSITION LAYER 
2.1. Physical Conditions in X-ray Skin Near Magnetic Flares 

The two-phase model was put forward by Haardt & Maraschi 
(1991, 1993) to explain the spectra of Seyfert Galaxies. Haardt 
et al. (1994) pointed out that observations are inconsistent with 
a uniform corona and introduced a patchy corona, where each 
“patch” is a magnetic flare (also referred as an active region). 
The key assumptions of the model are (1) during the flare, the 
X-ray flux from the active region greatly exceeds the disk in¬ 
trinsic flux, and ( 2 ) the compactness parameter I of the active 
region is large, so that the dominant radiation mechanism is 
Comptonization of the disk thermal radiation. We wish to ap¬ 
ply the same model to the X-ray spectra of GBHCs, and we 
employ the same assumptions. 

As discussed below, we argue for the existence of a transi¬ 
tion layer in the vicinity of an active coronal region (see Fig¬ 
ure 1). Since the flux of ionizing radiation is proportional to 
1/d^ X cos i oc d~^, where i is the angle between the normal of 
the disk and the direction of the incident radiation field and d is 
the distance between the active region and the position on the 
disk, the ionization state of the disk surface will vary across the 
disk. Consequently, only the gas near the active regions (with a 
radial size ~ a few times the size of the active region, situated 
directly below it) may be highly ionized. Most reprocessing 
of coronal radiation will take place in the transition regions; 
in addition, most radiation emitted by the disk that propagates 
through the active regions will have been emitted in the vicinity 
of the transition regions. Therefore, in this paper, we will only 
consider a one-zone model, consisting of the active region, the 
transition layer, and the underlying cold disk. 

The structure of the transition layer, i.e., its temperature, den¬ 
sity and ionization state are determined by solving the energy, 
ionization and pressure balance conditions. The first two condi¬ 
tions have been extensively treated in the lit erat ure, and we will 
follow these standard methods here (see § [2.2[ ). The pressure 
equilibrium condition is typically replaced by the constant gas 
density assumption (e.g., RFB, Zycki et al. 1994 and references 
therein). For the problem at hand, however, the equilibrium 
state of the transition layer is sensitive to the pressure balance, 
and thus we will attempt to take it into account. Accounting 
for the pressure balance leads to results that are substantially 
different from previous work, as discussed below. 

The X-ray radiation pressure on the transition layer is equal 
to F^/c, where is the flux produced by the magnetic flare. 


The compactness parameter of the active region, I, is defined as 


F^ctt^R 

^ ’ 

UleC^ 


( 1 ) 


and is expected to be larger than unity (e.g., Poutanen & Svens- 
son 1996, Poutanen, Svensson & Stern 1997). The size of the 
active region AR is thought to be of the order of the accretion 
disk height scale FI (e.g., Galeev et al. 1979). Here, H/Ris es¬ 
timated from the gas pressure dominated solution of Svensson 
& Zdziarski (1994; SZ94 hereafter). 


^ = 7.5 X 
K 

( 2 ) 

where a is the viscosity parameter. Mi = M/10Mq is the 
mass of the black hole, / is the fraction of accretion power 
dissipated into the corona (averaged over the whole disk), r is 
the radius relative to the Schwarzschild radius (r = R/Rg), 

J{r) = 1 — (3/r)^/^. For the case of the hard state of Cyg 
X-1, most of the bolometric luminosity is in the hard X-ray 
band (e.g., Gierlinski et al. 1997). Thus, most of the accretion 
energy must be dissipated directly in the corona, i.e., / ^ 1 
(Haardt and Maraschi 1991; Stern et al. 1995). The dimension¬ 
less accretion rate rh = rjM(? /L-^dd for Cyg X-1 seems to be 
around 0.05. Here, M is the accretion rate, p = 0.056 is the ef¬ 
ficiency for the standard Shakura-Sunyaev disk, and LEdd is the 
Eddington luminosity. Note that this definition of m is differ¬ 
ent by factor 77 from that used by SZ94 (i.e., rh ~ 17 x rhszgd)- 
Finally, for r = 6 , the X-ray flux is 

F^ - 4xl023/ai/i°M-®/'° ( ^ j (l-/)-i/i0ergcm-2sec-i 

( 3 ) 

To check whether assumption (1) of the patchy two-phase 
model is consistent for Cyg X-1 parameters, we estimate the 
intrinsic flux of the disk. 


Fd = 1.0 X lO^^M;^ ^ ( 0 ^) ~ ^ 

It is seen that the X-ray flux is indeed much larger than the 
intrinsic disk emission if 1 — / -C 1 and the compactness pa¬ 
rameter I 3 > 0 . 01 . 

The pressure of the disk surface layer before the occurrence 
of a flare (or, equivalently, far enough from the flare), assum¬ 
ing that the upper layer of the disk with Thomson optical depth 
Tx ^ few is in the hydrostatic equilibrium, is 

Po 

XTx j (1 - erg cm-3, (5) 

where r = 6 (SZ94). Near an active magnetic flare, the ratio 
of the incident radiation pressure to the unperturbed accretion 
disk atmosphere pressure is 

( 6 ) 

Since the radiation pressure from the active region greatly ex¬ 
ceeds the unperturbed thermal pressure, the transition layer will 
contract until the gas pressure P 

P ~ Px/c . 


( 7 ) 
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Fig. 1.— The geometry of the active region (AR; magnetic flare) and the transition layer. Magnetic fields, containing AR and sup¬ 
plying it with energy are not shown. Transition region is defined as the upper layer of the disk with Thomson depth of ~ few, where 
the incident X-ray flux is substantially larger than the intrinsic disk flux 


A better treatment would solve for the gas opacity in the tran¬ 
sition layer and thus the radiation force self-consistently, rather 
than simply using the ram pressure F^jc (Nayakshin 1998). 
However, inequality (Q) turns out to be sufficient to prove the 
main point of this paper. We thus move on to solve the energy 
and ionization balance equations for the transition layer. 

2.2. The Thermal Instability 


A general condition for a thermal instability was discovered 
by Field (1965). He argued that a physical system is usually in 
pressure equilibrium with its surroundings. Thus, any pertur¬ 
bations of the temperature T and the density n of the system 
should occur at a constant pressure. The system is unstable 
when 


/ 9Anet 
V dT 


p 


< 0 , 


( 8 ) 


where the “cooling function,” Anet, is the difference between 
cooling and heating rates per unit volume, divided by the gas 
density n squared. 

In ionization balance studies, it turns out convenient to define 
two parameters. The first one is the “density ionization param¬ 
eter” equal to (Krolik, McKee & Tarter 1981) 
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The second one is the “pressure ionization parameter”, defined 
as 


^ _ .frad 


( 10 ) 


where P is the full isotropic pressure due to neutral atoms, ions, 
electrons and the trapped line radiation. (This definition of 2 is 
the one used in the ionization code XSTAR, and is more appro¬ 
priate for studies of the instabilities than the original definition 
of 2 given in Krolik et al. 1981, who used nkT instead of the 
full pressure P. For temperatures typical of GBHCs, however, 
we found that the trapped line radiation was always a small frac¬ 
tion, e.g., ~ few percent of the total pressure, and therefore the 
two definitions of 2 are nearly identical). Krolik et al. (1981) 
showed that the instability criterion (|^ is equivalent to 


dT J 


< 0 


( 11 ) 


^net —0 


where the derivative is taken with the condition Anet = 0 sat¬ 
isfied, i.e., when the energy balance is imposed. In this form. 


the instability can be easily seen when one plots temperature T 
versus 2. 

We now apply the X-ray ionization code XSTAR (Kallman 
& McCray 1982, Kallman & Krolik 1986), to the problem of 
the transition layer. A truly self-consistent treatment would in¬ 
volve solving radiation transfer in the optically thick transition 
layer, and, in addition, finding the distribution of the gas den¬ 
sity in the transition layer that would satisfy pressure balance. 
Since radiation force acting on the gas depends on the opacity 
of the gas, this is a difficult non-linear problem. Thus, we defer 
such a detailed study to future work, and simply solve (using 
XSTAR) the local energy and ionization balance for an opti¬ 
cally thin layer of gas in the transition region. We assume that 
the ionizing spectrum consists of the incident X-ray power law 
with the energy spectral index typical of GBHCs in the hard 
state, i.e., F = 1.5 — 1.75, exponentially cutoff at 100 keV, and 
the blackbody spectrum from the cold disk below the transition 
layer. If the energy and ionization balance is found to be un¬ 
stable for this setup, the transition layer will also be unstable, 
because the instability is local in character. 

Note that it is not possible for the transition region to have 
a temperature lower than the effective temperature of the X- 
radiation, i.e., Tmin = {Fx/aY^^. The reason why the sim¬ 
ulations may give temperatures lower than Tmin for low val¬ 
ues of ^ is that in this parameter range XSTAR neglects certain 
non-radiative de-excitation processes, which leads to an over¬ 
estimate of the cooling rate (Zycki et al. 1994; see their sec¬ 
tion 2.3). Since we are using a one-zone model for the tran¬ 
sition layer, we use an attenuated X-ray flux (Fx) which rep¬ 
resents the surface-averaged value as seen by the transition re¬ 
gion, (Fx) = Fx/6 = 0.1Fx/&i, where & is a dimensionless 
number of order 10, bi = &/10, and Fx is the X-ray flux leav¬ 
ing the active region (see figure |^. Using equation (|]), 


- 5.0x10® 


0.05 


- 1/20 


^1/40 ^-9/40 j 


( 12 ) 

Figure ^ shows results of our calculations for several dif¬ 
ferent X-ray ionizing spectra. A stable solution for the tran¬ 
sition layer structure will have a positive slope, and also sat 


isfy the pressure equilibrium condition. As discussed in ^2.1 


P < Fx/c (i.e., > 1). In addition, if the gas is completely 

ionized, the absorption opacity is negligible compared to the 
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Fig. 2.— Gas temperature versus pressure ionization parameter S - the ionization equilibrium curves for parameters appropriate 
for GBHCs. The incident spectrum is approximated by a power law of photon index F, exponentially cutoff at 100 keV, and the 
reflected blackbody with equal flux and temperature T„iin (equation Values of the parameters are: F = 1.5, 1.75, 1.75, 1.7 and 
fcTmin = 200, 100, 200, 400 eV, corresponding to the fine solid, thick solid, dotted and dash-dotted curves, respectively. Because the 
ionization equilibrium is unstable when the curve has a negative slope, and because there exists no solution below Tmin (see text), 
the only stable solution is the uppermost branch of the curve with T ^ 10^ K. 


Thomson opacity. Because all the incident X-ray flux is eventu¬ 
ally reflected, the net flux is zero, and so the net radiation force 
is zero (note that the momentum of the incident radiation is re¬ 
flected deep inside the disk, so the radiation does apply a ram 
pressure equal to 2Fx/c to the whole disk, but not to the transi¬ 
tion layer). In that case the pressure P adjusts to the value ap¬ 
propriate for the accretion disk atmosphere in the absence of the 
ionizing flux Gee also Sincell & Krolik 1996), which is given 
by equation (^. Therefore, the pressure ionization parameter 
should be in the range 

1 < S < 1 X 10^; (1 - (13) 

With respect to the ionization equilibria shown in Figure (H), 
the gas is almost completely ionized on the upper stable 
branch of the solution (i.e., the one with T ^ 10^ K), and 
thus the pressure equilibrium for such temperatures requires 
S ~ Fx/cPq 3> 1, which is allowed according to equation 

In addition to the stable Compton equilibrium state, there is 
a small range in parameter space, with the temperature between 
100 and 200 eV, in which a stable thermal state exists . The 
presence of this region is explained by a rapid decrease in heat¬ 
ing (with T increasing), rather than an increase in cooling (cf. 
equation H). The X-ray heating rate decreases in the tempera¬ 
ture range 100 — 200 eV with increasing T since higher tem¬ 
peratures lead to higher ionization rates of ion species with ion¬ 
ization energies ~ kT. This larger degree of ionization reduces 
the X-ray opacity, and thus the heating rate as well. Note that it 
is highly unlikely that the transition region will stabilize within 
the temperature range 100 - 200 eV because the effective tem¬ 
perature, Tinin, is most likely above this temperature range. 

Due to the above considerations, it is very likely that the tran¬ 


sition layer is highly ionized in GBHCs in the hard state for 
Tx ^ 1. The upper limit of Tx can only be found by a more 
exact treatment. In addition, the transition layer may be heated 
by the same process that heats the corona above it, albeit with 
a smaller heating rate. Furthermore, Maciolek-Niedzwiecki, 
Krolik & Zdziarski (1997) have recently shown that the ther¬ 
mal conduction of energy from the corona to the disk below 
may become important for low coronal compactness parame¬ 
ters and substantially contribute to the heating rate of the tran¬ 
sition layer. Thus, the transition layer may be even hotter than 
found by photoionization calculations. 

Eventually, the X-rays are down-scattered and the radiation 
spectrum becomes softer as one descends from the top of the 
transition layer to its bottom. We can qualitatively test the gas 
ionization stability properties by allowing the ionizing spec¬ 
trum to be softer than the observed spectrum of GBHCs in the 
hard state. In Figure || we show two examples of such calcu¬ 
lations. The slope of the ionization equilibrium curve becomes 
positive everywhere above kT ~ 100 eV, so that these equi¬ 
libria are stable, and thus the gas temperature may saturate at 
T ~ Tmin, far below few keV, the appropriate temperature for 
the uppermost layer of the transition region. Thus, we know 
(see also that the transition layer should terminate at some 
value of Tx ~ few. 

We point out the similarity of our results to those of Kro¬ 
lik, McKee & Tarter (1981) for the emission line regions in 
quasars. These authors solved the ionization and energy bal¬ 
ance equations for optically thin clouds illuminated by a broad 
band quasar spectrum, and showed that the thermal instabil¬ 
ity exists if the pressure ionization parameter is close to unity. 
They found that there are two stable states for the line emitting 
clouds, one cold and one hot. The cold state corresponds to the 
gas temperature T ^ few x 10"^ K, where collisional cooling 
balances ionization heating. The hot state corresponds to tern- 
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Fig. 3.— Same as Figure ^ The thick solid curve is same as that in Figure^ and is appropriate for the hard state of a GBHC, whereas 
two other curves are relevant to the soft state in GBHCs, or at large depth in the transition layer (see Discussion). Values of the 
parameters are: F = 2.1, 2.1 and /cTmin = 200, 400 eV for the dotted and dashed curves, respectively. 


peratures T ^ 10^ K, when ionization heating decreases due 
to almost complete ionization of the elements and the Compton 
heating decreases since the gas is close to the Compton temper¬ 
ature Tcomp- Here, we found a similar instability for the X-ray 
skin near an active magnetic flare in accretion disks of GBHCs. 
The lower temperature equilibrium state is not allowed, how¬ 
ever, since the gas density and the X-ray ionizing flux in GB¬ 
HCs is larger than these quantities in quasar emission line re¬ 
gions by some ~ 10 — 12 orders of magnitude. Finally, we note 
that the thermal instability is not apparent in studies where the 
gas density is fixed to a constant value. Following Field (1965), 
we argue that the assumption of the constant gas density is not 
justified for real physical systems, and that one always should 
use the pressure equilibrium arguments to determine the actual 
gas density and the stability properties of the system. 

3. GLOBAL ADC MODELS WITH A TRANSITION LAYER 

3.1. Physical Setup 

We now explore how a transition layer affects the physical 
properties of the coronal gas and the spectrum of escaping radi¬ 
ation. Since a transition layer may occur for both a global ADC 
model and for a patchy coronal model, both models should 
be explored in detail. In this paper, however, we only study 
the global ADC model, and defer a self-consistent treatment of 
the patchy corona to a future paper. We expect, however, that 
the systematic trends found for the global ADC models should 
also occur for the patchy model, with the main difference being 
the maximum coronal temperature of the patchy model being 
higher. 

A self consistent treatment of the transition layer would solve 
for the ionization layer structure as a function of optical depth. 
However, we will make the simplifying assumption that the 
transition layer is completely ionized. The optical depth of the 


transition layer will be kept as a free parameter, and we will 
explore how sensitive the spectrum of escaping radiation is to 
Ttr. The two issues we will address are (1) for a given coro¬ 
nal optical depth (e.g., one thought to be appropriate for fitting 
the spectra of GBHCs), how does the maximum coronal tem¬ 
perature depend on the optical depth of the transition layer, and 
(2) how does the spectmm of escaping radiation, most impor¬ 
tantly the reprocessing features, vary with the optical depth of 
the transition layer. Of particular importance is the parameter 
range of Ttr such that the predicted spectrum is consistent with 
that of GBHCs, and whether these values of Ttr are consistent 
with the assumption of the transition layer being completely 
ionized. 

The model contains three regions: (1) A cold accretion disk, 
assumed here to have a temperature fcTbb = 150 keV, (2) the 
transition layer, situated directly above the cold disk, and (3) 
the corona, situated directly above the transition layer. Plane 
parallel geometry is assumed. 

We use the slab-geometry ADC model of Dove, Wilms, 
& Begelman (1997), which uses a non-linear Monte Carlo 
(NLMC) routine to solve the radiation transfer problem of the 
system. The free parameters of the model are the seed optical 
depth Te, (the optical depth of the corona excluding the contri¬ 
bution from electron-positron pairs), the blackbody temperature 
of the accretion disk and its compactness parameter, Zbb, and the 
heating rate (i.e., the compactness parameter), Ic, of the ADC. 
The temperature structure of the corona is determined numer¬ 
ically by balancing Compton cooling with heating, where the 
heating rate is assumed to be uniformly distributed. The 
pair opacity is given by balancing photon-photon pair produc¬ 
tion with annihilation. Reprocessing of coronal radiation in the 
cold accretion disc is also treated numerically. For a more thor¬ 
ough discussion of the NLMC routine, see Dove, Wilms, & 
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Fig. 4.— The maximum temperature of the corona as a function of the optical depth of the transition layer. For all models, we 
assume that the intrinsic compactness parameter of the disk is Zbb = 0.01 and the disk temperature is fcTbb = 150 eV. For all models, 
the maximum temperature is reached for Ic ^ 2 — 5. 
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Fig. 5.— The predicted spectrum for various values of the transition layer optical depth. From top to bottom, Ttr = 10, 5, 2.5, and 
1 . 0 . 
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Begelman (1997). The transition layer is treated identically to 
the corona, accept here the heating rate is set to zero. Therefore, 
the transition layer, numerically modeled using 8 shells, each 
with equal optical depth dr = Tt^/S, will obtain the Comp- 
tonization temperature due to the radiation field from both the 
corona and the accretion disk. 

3.2. Maximum Coronal Temperature 

As discussed by Dove, Wilms, & Begelman (1997), for a 
given total optical depth, there exists a maximum coronal tem¬ 
perature, above which no self-consistent solution exists. Rais¬ 
ing the compactness parameter of the corona to a value higher 
than that corresponding to the maximum temperature gives rise 
to a higher optical depth (due to pair production), causing more 
reprocessing, subsequent Compton cooling, and ultimately a 
lower coronal temperature. Here, as much as 80-90% of the 
incident X-ray flux is re-radiated at the disk temperature, i.e., 
the X-ray integrated albedo is 0.1-0.2. A completely ionized 
transition layer will increase the albedo of the cold disk. Ac¬ 
cordingly, the amount of energy in the reprocessed blackbody 
(or Comptonized blackbody) should become smaller with in¬ 
creasing Thomson optical depth of the transition layer, and 
this Compton cooling via reprocessed radiation should be less 
efficient and therefore higher maximum coronal temperatures 
(as compared to models with no transition layer) should be al¬ 
lowed. 

In Figure we show how the average coronal temperature 
varies with the optical depth of the transition layer, Ttr. In Fig¬ 
ure ^ we also show the resulting broad band spectra for the 
model parameters tested. For all models, we chose Zbb = 0.01 
in order to be consistent with the definition of the disk com¬ 
pactness parameter, i.e., Zbb = /{meC^), according to 

which 

/bb = 0.03(^^j (14) 

Further, other parameter values are Ic = 2, kTbi, = 150 keV, 
and Tc = 0.3. These parameters correspond to the model pro¬ 
ducing the maximum coronal temperature. In contrast to mod¬ 
els in which Ttr = 0, the coronal temperature for a given value 
of Tc is not simply a function of Ic/lhb- To see this, consider 
the case where Ttr 1- Here, the albedo of the disk is essen¬ 
tially unity, and therefore all of the soft photons emitted will 
be from the intrinsic flux of the disk (no reprocessing). There¬ 
fore, setting ^bb ^ 1 yields the maximum coronal tempera¬ 
tures possible. Note that the maximum temperature levels out 
as Ttr ^ c)o. Although, in this limit, there is no reprocessing 
of hard X-rays in the cold disk, there is still “reprocessing” in 
the transition layer. As Ttr increases, more coronal radiation 
is down-scattered to the Compton temperature of the transition 
layer, which is kTtr ~ 1 — 4 keV. Even at these temperatures, 
Compton cooling of this “reprocessed” radiation in the corona 
is very efficient. 

It is interesting to note that, only for Ttr ^10, the coronal 
temperature is high enough such that the corresponding spec¬ 
trum of escaping radiation is hard enough to describe the ob¬ 
servations of Cyg X-1. (The canonical value of the photon 
power law of Cyg X-1 is F = 1.7; for Tc = 0.3, this power 
law corresponds to kTc ~ 150 keV). It is probably unphysi¬ 
cal, however, to assume the transition layer is completely ion¬ 
ized for such large optical depths. In fact, the numerical model 
for Ttr = 10 predicts a temperature of fcTtr 500 eV near 
the bottom of the layer. Furthermore, as is seen in Figure ^ 


the radiation spectra do become substantially softer deeper in 
the transition layer, and, using our simple experimentation with 
softer ionizing spectra, shown in Figure we can expect that 
the gas will become stable and may be at a temperature lower 
than the Compton temperature for Ttr as small as 2 — 3. There¬ 
fore, even with the advent of transition layers, it still appears 
unlikely that a global slab geometry ADC model can have self- 
consistent temperatures high enough to reproduce the observed 
hard spectra of Cyg X-1 and other similar BHCs. This does not 
rule out ADC models in which the corona is patchy, e.g., con¬ 
taining several localized active regions such as magnetic flares. 
The patchy geometry leads to higher coronal temperatures due 
to less reprocessed soft flux re-entering active regions (Pouta- 
nen & Svensson 1996), so that lower Ttr may be sufficient to 
explain the Cyg X-1 spectrum. 

3.3. Reprocessing Features in the Spectrum of Escaping 
Radiation 

As shown in Figure ||, the reprocessing features in the spec¬ 
trum of escaping radiation depend sensitively on the optical 
depth of the transition layer. In Figure ^ we show how the 
internal spectrum (averaged over all angles) varies throughout 
the transition layer and corona. Note that most of the hard X- 
rays do not penetrate through the transition layer, and that the 
spectrum gets softer as it approaches the cold disk. On the other 
hand, the thermal radiation emitted by the disk and the 10-50 
keV Compton reflection hump are broadened by Compton scat¬ 
terings as the disk radiation diffuses upward through the transi¬ 
tion layer and the corona. The Iron Ka line, small to start with 
due to the small amount of reprocessing of coronal radiation, 
is completely smeared out by the time the radiation escapes the 
system. No line photons are created in the transition layer itself, 
because we found that the Compton equilibrium state typically 
resides at the ionization parameter ^ 10'^, whereas no fluo¬ 

rescent iron line emission is produced for ^ > 5 X 10^ (Matt, 
Fabian & Ross 1993, 1996). 

Therefore, in agreement with Bottcher, Liang, and Smith 
(1998), we conclude that a highly ionized transition layer can 
reduce the reprocessing features in the spectrum of escaping 
radiation emanating from slab geometry ADC models. Previ¬ 
ous work has used the lack of observed reprocessing features 
to argue against the applicability of slab geometry ADC mod¬ 
els, and argued for a disk/corona configuration where the disk 
has a small solid angle relative to the coronal region, such as a 
“sphere-fdisk model (Dove et al. 1998, Gierlinski et al. 1997, 
Poutanen, Krolik & Ryde 1997). A slab geometry ADC model 
containing a transition layer, however, is another possibility. 
Although, as discussed above, these models could not obtain 
high enough temperatures to be able to predict spectra as hard 
as the observed spectra of GBHCs unless Ttr ^ 10, a model 
with a patchy corona and underlying transition layers appears 
likely to be able to predict hard enough spectra with a relatively 
small amount of reprocessing features. In a forthcoming paper, 
we plan to investigate such a model. 


4. DISCUSSION 

We have shown that the transition layer in the vicinity of tran¬ 
sient flares for ADC models of GBHCs must be highly ionized 
and very hot. Specifically, due to a thermal instability, the only 
stable temperature of the transition layer is the local Compton 
temperature {kT ~ few keV). In fact, even for global ADC 
models of GBHCs, such a transition layer is found to be likely. 
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Fig. 6.— Internal spectrum for several shells within the transition disk and corona. From bottom to top, Ttr(z) = 1.25, 2.5, 3.75, and 
5.0, where Ttr{z) is the Thomson optical depth from the cold-disk/transition layer interface {z = 0) to a height z. The uppermost 
spectrum is the internal spectrum within the corona. For this model, rtr = 5.0 and Tc = 0.3. 


However, we have not accurately modeled the full radiative 
transfer problem through the transition layer at this time, and 
therefore could not determine its vertical optical depth, which 
was treated as a free parameter. 

Due to the transition layer, a larger fraction of incident X- 
rays are Compton reflected back into the corona without being 
reprocessed by the cold disk. In addition, for Ttrans ^ 1, the 
predicted reprocessing features as well as the thermal excess 
should be substantially smaller (as is found in this paper for 
the global ADC model) than that of previous ADC models in 
which the transition layer was not considered. This reduction 
of the reprocessing features is crucial for the model being con¬ 
sistent with the observations of GBHCs (e.g., Gierlinski et al. 
1997, Dove et al. 1998). 

This reduction in reprocessing yields a lower Compton cool¬ 
ing rate within the corona, and higher coronal temperatures than 
previous ADC models are allowed. For global ADC models 
with Tc ~ 0.3, we find that the coronal temperature can be as 
high as ^ 150 keV if the optical depth of the transition layer is 
Ar ^ 10. However, a completely ionized transition layer with 
Ttr ^ 10 is most likely physically inconsistent since the bot¬ 
tom layer was found to be too cold to realistically stay highly 
ionized. Therefore, global slab-geometry ADC models are still 
problematic in explaining the observations of GBHCS. Nev¬ 
ertheless, these global ADC results are very encouraging, and 
a slab geometry ADC model containing a patchy corona (e.g., 
individual ARs) with underlying transition layers should be rig¬ 
orously studied. Here, due to the lower amount of reprocessed 
radiation within the ARs as compared to the global model, mod¬ 
els with coronae hot enough to reproduce the observed spectra 


of GBHCs may be allowed for more reasonable transition layer 
optical depths. 

We have considered only radiative cooling mechanisms for 
the transition layer in this paper. It is possible that a wind 
is induced by the X-ray heating. However, the maximum gas 
temperature obtained due to the X-ray heating is the Comp¬ 
ton temperature (^ 10® K). Therefore, as shown by Begelman, 
McKee & Shields (1983), a large scale outflow cannot occur 
for R ^ 10''^i?g. On the other hand, a local uprising of the 
gas is still possible. The maximum energy flux due to this pro¬ 
cess is Fev ~ Pcs, where Cg is the sound speed in the transi¬ 
tion region. Since P ^ Fx/c, and Cs ^ 3 x 10“®c, we have 
Fev/Fx ^ 3x10“®. Therefore, a wind or any of/zer mec/zanica/ 
process cannot cool the gas efficiently, and thus it is justifiable 
to consider cooling via emission of radiation only. 

In a forthcoming paper, we will discuss the implications of 
the instability for the case of AGN. As shown in Figures ^ and ^ 
there is a stable ionization equilibrium state below T ~ 3 x 10® 
Kelvin. Furthermore, due to increasing opacity to the soft disk 
radiation, we find that equilibria below T ~ 10® K are not pos¬ 
sible. We found that the existence of this stable region, and 
the unstable region above 3 x 10® K, is very much independent 
of the details of the incident spectra (indeed. Figures g and || 
were not even intended for the AGN case). Thus, it is possible 
that the differences in the intrinsic X-ray spectra, the ionization 
state of the reflector and the strength of the iron line in AGNs 
and GBHCS can be explained by the different end points of the 
thermal ionization instability described here. 
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